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We study electronic properties of a superconducting topological insulator whose parent material 
is a topological insulator. We calculate the temperature dependence of the specific heat and spin 
susceptibility for four promising superconducting pairings proposed by L. Fu and E. Berg (Phys. 
Rev. Lett. 105, 097001). Since the line shapes of temperature dependence of specific heat are 
almost identical among three of the four pairings, it is difficult to identify them simply from the 
specific heat. On the other hand, we obtain wide varieties of the temperature dependence of spin 
susceptibility for each pairing reflecting the spin structure of Cooper pair. We propose that the 
pairing symmetry of superconducting topological insulator can be determined from measurement of 
Knight shift by changing the direction of applied magnetic field. 



I. INTRODUCTION 

Topological insulators (TIs) arc new state of mat- 
ter supporting masslcss Dirac fcrmions on the surface, 
which are characterized by non-zero topological num- 
bers defined in the bulkii^ Due to the surface Dirac 
fcrmions, TIs have potential to exhibit rich transport 
and electromagnetic response properties, which can be 
applicable for future devices. The superconducting ana- 
log of TI are topological superconductors pT— which have 
Majorana fcrmions 7 - on the surface as Andreev bound 
states (ABSs). In these materials, topological invari- 
ants can be defined in the bulk Hamiltonian. There are 
several types of topological superconductors, e.g., chi- 
ral p-wave superconducting state in Sr2Ru04,— ~— helical 
superconducting state realized in non-centrosymmetric 
superconductors! 12 i 13 The realization of a topological su- 
perconductor is of particular interest in the viewpoint of 
quantum devices and quantum computations.— ~— 

Recently, carrier-doped TI Cu x Bi 2 Se3 is revealed to 
be a superconductor^ Hereafter, we refer to a super- 
conductor based on TI as a superconducting topolog- 
ical insulator (STI). In the tunneling spectroscopy,— 
Cu,;Bi2Sc3 shows a zero-bias conductance peak (ZBCP). 
This means that Cu T Bi 2 Sc3 can be regarded as a topo- 
logical superconductor since the ZBCP signifies the ex- 
istence of gaplcss ABSs^i~— on the surface, which is a 
direct consequence of topological superconductivity. In- 
terestingly, it has been clarified that the STI supports 
anomalous ABSs different from those of other topolog- 
ical superconductors , 29 ' 37 ! 38 and the resulting transport 
property also becomes anomalous! 29 ' 38 In this sense, STIs 
are new type of topological superconductor, and has at- 
tracted much interest. Moreover, there are several ex- 
perimental results supporting a generation of STI by the 
proximity effecti 39 ' 40 

Nowadays, there are a lot of relevant studies about 
Cu E Bi2Sc3 . 32 ' 33 i 41 ~ — However, up to now, the symme- 
try of superconductivity of Cu x Bi2Sc3 still remains un- 
known. Although the specific heat has been measured, 
it is difficult to establish the superconducting symmetry 
only from the data of specific heat. More careful analysis 
with the help of the microscopic calculations is needed. 



In order to clarify the superconducting symmetry, it is 
useful to analyze spin susceptibility besides specific heat 
since the spin susceptibility is directly related to spin 
structure of superconducting pairing. 

In this paper, we clarify the temperature dependences 
of specific heat and spin susceptibility for the possible 
superconducting pairings. Based on these quantities, it 
is possible to determine the pairing symmetry in a STI. 

The paper is organized as follows. In Sec. UH we show 
the model Hamiltonian of STI, and give the energy spec- 
tra for the possible superconducting parings. The nu- 
merical results and discussion on the temperature de- 
pendences of the specific heat and spin susceptibility are 
shown in Sec. IIII1 and Sec. IIV[ respectively. We compare 
our results of specific heat with the experimental data.— 
In Sec. El we summarize our results and propose how to 
experimentally determine the superconducting symmetry 
of STI. 



II. MODEL 

For a model of STI, we start with the BdG Hamiltonian 
proposed in Refliil. 



H(k) = H Q (k)r z + A e r x , 



(1) 



where t = 1, 2, 3, 4 represent the type of the pair poten- 
tial. The normal part of the Hamiltonian Ho(k) is low- 
energy effective model of a topologocal insulator, based 
on the k ■ p theory given by 



Ho{k) 

m(k) 
c(fe) 



c(fc) + m(k)a x 

+v z k z a y + v(k x s y - k y s x )a z , (2) 
m a + mifc 2 + m 2 (kl + fc 2 ), (3) 



= -Li + Cl kl+c 2 {k 2 x + k 2 y ). 



(4) 



Oi and Tj (i = x, y, z) are the Pauli matrices in spin, 
orbital and Nambu spaces, respectively. The basis of the 
orbitals are effective p z orbitals constituted from p z or- 
bitals of Se and Bi p z orbitals on the upper and lower side 
of the quintuple layer, as shown in Fig. [TJ Hereafter, we 
call this basis "orbital basis" . On the other hand, we refer 
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to the basis diagonalizing Ho(k) as "band basis", which 
is introduced in Sec. IIVBI In this model, the normal 
part H (k) is equivalent to the model proposed in Refill 
and Refl50l under the unitary transformation. In the fol- 
lowing, we use the tight-binding model which is equiva- 
lent to the above model at low energy. We consider the 
hexagonal lattice where two-dimensional triangular lat- 
tices stack along the c-axis i 33 ' 37 Then, the tight-binding 
Hamiltonian is obtained by the following substitution in 
the k ■ p Hamiltonian Eqs. ©-flU- 



2 . V3k x a k y a 

sin — - — cos -f- = f x (k), 



2 



(5) 



2 , y/3k x a . 

— (cos — - — sin + smsjoj = J y (k), 



k 2 



■ smk z c = f z (k), 



(1 - cosfc 2 c) = f±(k), 



= /||(fc), 



\/3k x a k v a 
2 cos — - — cos — | cos k y a) 



(6) 
(7) 

(8) 
(9) 



where a and c are the lattice constants. In this hexagonal 
lattice, the primitive lattice vectors are (\/3o/2, a/2, 0), 
(0,a, 0), (0,0, c), though the actual crystal structure is 
not hexagonal but rhombohcdral onei 49 i 50 This simplifi- 
cation does not affect the low energy excitations. Then, 
the normal part of the Hamiltonian is summarized as fol- 
lows. 



H (k) 


= c(k) + m(k)a x 








+ (a x (k)s y - a y 


(k)s s )<r z + b(k)a y , 


(10) 


c(fe) 


= -fj, + ci/j_(fe) H 


-c 2 /||(fe), 


(11) 


m(k) 


= m + mi/j_(fe) 


+ m 2 f\\(k), 


(12) 


a x,y(k) 


= vf x>y (k), 




(13) 


b(k) 


= v z f z (k). 




(14) 



Here, we choose the chemical potential fi = 0.5 eV, since 
the chemical potential measured from the surface Dirac 
point is 0.4^0.5 cV according to Reffill. We use the val- 
ues of the parameters, C2, mo, 777,2, v as given in RefUfJ. 
On the other hand, for cj, m 1 , v z , we choose the dif- 
ferent values given in Ref|50j, which involve the hopping 
along the c-axis. Since the parameterization performed 
in ReffHol is based on the dispersion at around T-point, 
the difference of the dispersion near the zone boundary 
between the first principle calculation in Refisol and our 
tight-binding model is considerably large. However, the 
Fermi surface becomes a cylindrical one if we use the 
same parameters given in RefUFJ although the correct 
shape of the Fermi surface is an oval spherical one. Thus, 
we choose the values of ci, mi, v as ci/c 2 = 0.024, 
toi/c 2 = 0.216, v z /c — 0.32 (eV) to fit the energy disper- 
sion for the T-Z direction obtained in ReflHol. These pa- 
rameters give the oval spherical Fermi surface consistent 



with the first principle calculation. This parameteriza- 
tion is crucial since specific heat and spin susceptibility 
in actual Cua;Bi2Se3 can not be reproduced if we use a 
cylindrical Fermi surface. In addition, to obtain the topo- 
logical superconductivity in three dimensions, the correct 
Fermi surface topology is needed 

Next, we consider the pair potential. We assume that 
the pair potential is independent of momentum since the 
present material is not a strongly-correlated system4£ 
In this case, these pair potentials are classified into four 
kinds of irreducible representation for D%d point group. 
The matrix forms of the pairings Ai, A2, A3 and A4 are 
shown in the first column of Table. |U Ai and A3 are 
spin-singlet intra-orbital pairings, whereas A2 and A4 
are spin-triplet inter-orbit pairings in the orbital basis. 

We diagonalize the BdG Hamiltonian Eq.(fT]). We ob- 
tain four branches of the bulk spectrum E 1 (7 = 1, 2, 3, 4) 
for each pairing, 



Ex{k) = V£ 2 (fc) + 2^ 2 (fe)c 2 (fe) + C 2 (fc)A 2 , (15) 



E 2 (k) = Je(k) - 2^if{k)c 2 {k) + C 2 (fc)A 2 , (16) 



E 3 (k) = -\^e(k) + 2vV(fc)c 2 (fc) + C 2 (fc)A 2 ,(17) 



Ei{k) 



e(k) - 2V?7 2 (fe)c 2 (fc)+C 2 (fc)A 2 ,(18) 



with 



V 2 (k) 

e(k) 



m 2 {k) 
V 2 (k) 



Va 2 x {k) 
c 2 {k) 4 



ha 2 (fc) + 6 2 (fc), 
A 2 . 



(19) 
(20) 



The difference of energy gap structure in each pairing 
comes from ( 2 (k), 

(21) 
(22) 
(23) 
(24) 



Ax 


C 2 (k) 


= 0, 


A 2 


C 2 (k) 


= m 2 (k) 


A 3 


C 2 (k) 


= m 2 {k) 


A 4 


C 2 (k) 


= m 2 {k) 



b 2 (k), 

a 2 Jk). 



Energy gap structure of Ai is isotropic full gap, which 
is the same as the conventional BCS superconductors. 
In other cases, because of the presence of C, energy gap 
is modified from BCS gap structure. A2 is anisotropic 
full gap pairing. In the case of A3 and A4, energy gap 
has point nodes. Point nodes for A3 are on the poles. 
In the case of A4, point nodes appear on the fc H -axis. 
Although, in general, A4 is a linear combination of Aa y s x 
and /S.<j y Sy, we can choose A4 = a y s x without losing a 
generality. Energy gap of E 7 is influenced by the spin- 
orbit interaction v. To elucidate the role of the spin orbit 
interaction, we also consider the case of v = 0. In this 
case, E 1 for Ai has a full gap, E^s for A2 and A4 have 
line nodes on equator, and E 7 for A3 is gapless. 



III. SPECIFIC HEAT 

In this section, we calculate the specific heat below T c 
for each pairing symmetry. The specific heat is given by 
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FIG. 1: (color online) Two p^-orbitals in the quintuple layer 
of Bi 2 Se 3 . 
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TABLE I: The irreducible representation, the spin state, the 
orbital state and the energy gap structure in each pairing 
symmetry. In the bracket we denote the gap structure for 
v = 0. 



two different values. One is a = ascs, and the other 
is a = a where C S {T)/T at T = 0.53T C = T becomes 
equal to that for the normal state C n (T)/T ~ C n {T c )/T c 
as observed in the specific heat measurement^ 



A. Isotropic full gap Ai 




0.0 0.5 1.0 

T/T c 

FIG. 2: (color online) Temperature dependence of specific 
heat for Ai with a = 1.94(blue solid line) and a = 1.76(green 
dash line) . Black solid line shows specific heat for nomal state. 
Red dotted circles show the experimental data by Refl43l 
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2/3 
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E 

k~, 

N ^ 



k 7 



dE 7 {k) 
dE 7 (k) 



E 2 ,{k) 



dE^ 



(k) 



d/3 



El{k) + pE^k) 



dA dE 7 (k) 
5/3 OA 



(25) 



where N is the number of unit cells and (3 is 1/fceT with 
the Boltzmann constant k B and temperature T. We as- 
sume the temperature dependence of the pairing poten- 
tial is the scaled BCS one, A(T) = (a/a B cs)A B cs(7 1 ). 
The magnitude of a gives the ratio of A(T = 0) to T c , 
i.e., a = A(T = 0)/(k B T c ). This model is dubbed as a- 
modelj^i For Abcs(T), w e use the following phenomeno- 
logical form^ 



(26) 



with aecs = 1-76. 

Since a is the material-dependent parameter and it 
often deviates from the BCS value a-e.es = 1-76, we use 



In Fig. [21 we show the temperature dependence of 
C s for a = 1.76 (green dash line) and a = cto = 1.94 
(blue solid line). In the case of Ai, the energy spectrum 

is given by E(k) = ±j£±(k) + A 2 , where £±{k) is the 



dispersion of normal state. Therefore, the energy gap 
structure becomes isotropic s-wave one. Thus, the spe- 
cific heat near T = shows an exponential behavior. If 
we choose a = cvbcSj Tq and the magnitude of specific- 
heat jump are smaller than those of the experiment. To 
fit the experimental data, we choose a = ao = 1.94. 
Then, to satisfy the entropy balance relation, 



dT Cs(T)-C n (T) = 
T 



(27) 



the magnitude of the specific-heat jump at a — 1.94 be- 
comes larger than that for a = a?BCS- 

In the case of a = ao, the magnitude of the specific 
heat jump and the line shape are similar to those of the 
experimental ones. We note that the analysis performed 
in Ref|H is based on isotropic s-wave gap and therefore 
the obtained value of ao is almost the same. On the other 
hand, the value of a can be also estimated from the upper 
and the lower critical field in ReflS The estimated value 
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is a = 2.3 = a c . Therefore, the value of a for Ai 
deviates from that of a c . However, if we add a small k- 
dependent term allowed in A lg representation to Ai, the 
magnitude of specific heat jump for a = aecs could be 
small and the values of ao becomes large, then, ao = 2.3 
might be obtained. 

B. Anisotropic full gap A 2 




T/T c 

FIG. 3: (color online) Temperature dependence of specific 
heat for A2 with a = 2.09. The dotted circles show the 
experimental data by Refl43l 



In Fig. [3l we show the temperature dependence of 
C S (T) for a = 1.76 (green dash line) and a = ao = 
2.09 (blue solid line). Since the energy gap structure 
is fully gapped, the exponential behavior appears near 
T = as in the case of Ai. On the other hand, the 
magnitude of specific heat jump for a = «bcs is smaller 
than that for Ai owing to the anisotropy of the energy 
gap. Therefore, to reproduce the experimental data, we 
need larger value of ao compared with the case of Ai, 
ao = 2.09. This value is closer to a c = 2.3 than that for 
Ai. The magnitude of the specific heat jump and the line 
shape for a = ao are similar to those of the experimental 
ones. 



C. Point nodes at polar A3 



of v = 0, the energy dispersion for A3 is given by E(k) — 
±Y / m 2 (fc) + b 2 (k)± y/c(k) 2 + A 2 . This energy spectrum 
becomes gapless when m 2 (k) + b 2 (k) = c(k) 2 + A 2 . Pa- 
rameters of STI satisfy the following relations. 

mo-// 2 -A 2 < 0, (28) 
m\-c\ > 0, (29) 
m\ - c\ > 0. (30) 

In this case, the energy spectrum becomes gapless near 
the Fermi surface in any direction of k. Thus, C S {T)/T 
with v — is T-indcpendent. In the presence of v, these 
gapless energy spectra still remain in the direction of k z - 
axis and point nodes are formed since a x , y {k) = in 
this direction. In the directions other than k z , the en- 
ergy gap is generated by the spin-orbit coupling v, but 
the gap is comparatively smaller than the other pairings. 
Therefore, the T-dependence of C S (T)/T for a = aecs 
remains relatively small. This is the reason why the spe- 
cific heat jump is small for A3 compared to other pair 
potentials. If we use a = ao = 2.74, we can adjust 
the specific heat jump similar to the experimental one 
and [C s (T )/T ]/[C n (T c )/T c ] becomes equal to the unity. 
However, the line shape does not reproduce the experi- 
mental data. In addition, the value of ao is much larger 
than the experimental value a c = 2.3. 




FIG. 4: (color online) Temperature dependence of specific 
heat for A3 with a — 2.74. The dotted circles show the 
experimental data by Refl43l 



In the case of A3, the energy dispersion has point 
nodes along fc z -axis. Thus, C S (T)/T has T 2 -bchavior 
near T = 0. The magnitude of the jump for a = aecs is 
the smallest among the four kinds of the pair potentials 
considered in this paper. This small jump is originated 
from the gapless nature of this pair potential. In the case 



D. Point nodes on equator A4 

In the case of A4, the energy spectrum has point nodes 
along fcy-axis. Therefore, C S (T)/T has T 2 -behavior near 
T = as in the case of A3. However, in the case of 
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A4, the energy spectrum does not become gapless even 
when v is absent. Thus, the coefficient of T 2 is smaller 
than the case of A3 for a = q;bcs , an d the magnitude of 
the specific heat jump is larger than that for A3. As a 
result, the line shape with a = = 2.42 is considerably 
close to the experimental one as compared with the case 
of A3. The obtained value of a — 2.42 is the most close 
to the experimental one, a c , among the four kinds of the 
pairing symmetry considered in this paper. 




FIG. 5: (color online) Temperature dependence of specific 
heat for A4 with a — 2.42. The dotted circles show the 
experimental data by Refl43l 



Here, we summarize the results of the specific heat. We 
have calculated the specific heat for a = cvbcs j <*o in each 
pairing symmetry. For a = ao, we obtain the line shapes 
similar to the experimental one in the cases of Ai, A2, 
and A4. The obtained values of ao are a = 1.94, 2.09, 
2.74, and 2.42 for Ai, A 2 , A 3 , and A 4 , respectively. The 
values of a for A2 and A4 are closer to the experimental 
one (ev c = 2.3), as compared with the other pair poten- 
tials. 



IV. SPIN SUSCEPTIBILITY 

From the temperature dependence of spin susceptibil- 
ity, one can determine the spin structure of Cooper pairs. 
Namely, for a spin-singlet superconductor, the spin sus- 
ceptibility along any direction decreases for T < T c and 
vanishes at T = 0. On the other hand, for spin-triplet su- 
perconductor, only the spin susceptibility parallel to the 
direction of d-vector decreases and vanishes at T = 0, and 
the spin susceptibility perpendicular to (i-vector is inde- 
pendent of T. However, the temperature dependence of 
spin susceptibility of a STI is not simple because spin- 
singlet and spin-triplet components can mix in the band 



pair potential 


effects of SOI 


Ax 


= A 


Van Vleck 


A2 


= A<JySz 


Van Vleck 
rotation of d-vector 


A 3 


= Act? 


Van Vleck 
induced spin-triplet 






Van Vleck 


A 4 


— Aa y s x 


rotation of d-vector 
induced spin-singlet 



TABLE II: Summary of spin susceptibility. Effects of spin- 
orbit interaction (SOI) is Van Vleck susceptibility, rotation 
of d-vector, and to induce spin-singlet and spin-triplet pair 
potentials. 



basis due to the spin-orbit interaction depending on the 
pair potential. 

Nevertheless, we show here that it is possible to deter- 
mine spin structure of STI, even if spin-orbit interaction 
is present. The temperature dependences of spin suscep- 
tibility with each possible pairing are different from each 
other. For Ai, the spin susceptibility along any direc- 
tions decreases as temperature decreases, since Ai is a 
spin-singlet pair potential in the band basis. On the other 
hand, that along z-axis for A3 is independent of temper- 
ature, although those along x- and y-axes decrease. Spin 
susceptibilities with A2 and A4 along the d- vector (d \\ z 
for A 2 and d \\ x for A 4 ) decrease. Those perpendicular 
to the (i-vector are almost independent of temperature. 
The temperature dependences of spin susceptibility are 
summarized in Tabic UTTl 

Now we comment on effects of spin-orbit interaction 
for spin susceptibility. There are three effects. The first 
one is Van Vleck susceptibility, which comes from inter- 
band (off-diagonal) matrix elements. Van Vleck suscep- 
tibility can appear in multiband systems with spin-orbit 
interaction. This leads to a non-zero value of spin sus- 
ceptibility at T = (See Appendix I A The second one 
is the rotation of (i-vector by the unitary transformation 
from the orbital basis to the band basis. Then, the d- 
vector in band basis is not parallel to the Zeeman mag- 
netic field, even when the (i-vector in the orbital basis is. 
This also induces a non-zero value of spin susceptibility 
at T = 0. Additionally, the spin susceptibility perpendic- 
ular to d- vector in the orbital basis also decreases slightly 
for T < T c . This behavior shows up in the case of A 2 
and A4. The third one is a generation of spin-singlet 
and spin-triplet pair potentials in the band basis from 
spin-triplet and spin-singlet ones in the orbital basis, re- 
spectively. In the following sections, we shall discuss the 
temperature dependence of the spin susceptibility in each 
pairing. 



A. Kubo formula for spin susceptibility 

First, we give the Zeeman term in a STI and Kubo 
formula for spin susceptibility. The Zeeman term Hz(k) 
is given by 



H z (k)= Y, 



E 

x,y,z fi—0,x,y,z 



(31) 



where fiB is Bohr magneton, hi is the i-th component 
of Zeeman field, gi^ is g-factor of the parent topological 
insulator, and (To is 2 x 2 identity matrix in the orbital 
space. The spin susceptibility along i-axis is given by the 
Kubo formula: 



Xi 



1 E f(E (k))-f( Ep ^ 



q)) 



iO 



{a\ Si \p){p\ 9 -f Si a»\a). 



(32) 



In actual calculation, we set g^ to be that of Bi2Se3 
9x0 = 9 y o = -8.92, g z0 = -21.3, g xx = g yx = QM,g zx = 
—29.5. The other g-factors are chosen to be zero. Tem- 
perature dependence of A is same as that used in the 
calculation of specific heat with a = 2.3. 



B. Spin structure in the band basis 

In the band basis, spin-singlet and spin-triplet pair po- 
tentials can mix each other because of spin-orbit interac- 
tion . Due to this, it is rather difficult to understand the 
temperature dependence of spin susceptibility. In order 
to clarify the spin structures of pair potentials, it is nec- 
essary to introduce the band basis where the normal part 
of Hamiltonian is diagonalized. First, we diagonalize the 
Hamiltonian Ho(k) of the normal state as 



H (k)uj(k) = e 7 (fe)it 7 (fe), 



(33) 



with band index 7. By using the unitary matrix U(k) 
given by U(k) = (ui(k), 112(h), ■ ■ ■ ), the pair potential 
A = ^ d^s^ is transformed as 



rfo(fc) and d(k) denote spin singlet-component of pair 
potential and <i-vector of spin-triplet component of pair 
potentials in the band basis, respectively. It is noted 
that s M does not change by this unitary transformation in 
inversion-symmetric systems. The corresponding Hamil- 
tonian is expressed as 



U\k)H(k)U(k) 



diag(ei(fe) - fJ,,e 2 {k) 
+ y^,d fl {k)s fl T x . 



(35) 




FIG. 6: (color online) Vector field plot of the conduction-band 
component of d(k) for A2 [(l,l)-component of Eq. (|37[) ]. 



In the following, we show the d-vectors in the band basis 
in the lowest order of k. The detailed derivation of d l _ l (k) 
is shown in Appendix. In the case of Ai, d^(k) is the 
same as in the orbital basis: do(k) = A,d(k) = 0. For 
other cases, we have 
A 2 : 



do(k) = 0, 

d(k) = A(^,^,^a z 
\m m |to | 



sgn(m )CT !/ 



As 



d (k) = Aa x , 
d(k) = Aa z 



vky vkx 



|to | ' \m \ 



,0 



A 4 



do(fc) 
d{k) 



Avky v z k z _ 
c a 

to too 

v z k z _ 



mo 



a y ,0, ■ 



vk x 
to 



(36) 
(37) 

(38) 
(39) 

(40) 
(41) 



Here, cr, is Pauli matrix denoting band index, i.e., a z = 
1 for the conduction band, a z = — 1 for the valence 
band. To illustrate the d(k) of Eq.([37|) in the conduc- 
tion (valence) band, we plot the (1, l)-component ((2, 2)- 
componcnt) of the d vector in Fig 15] (FigJT]). Those of 
Eqs. ((391) and (gTJ are also shown in Figs. [HJ O [lOl and 
II H respectively d^(k) is useful to understand the tem- 
perature dependence of xu as we will see in the following. 



C. Isotropic full-gap Ai: Van Vleck susceptibility 

Fig. [T3] shows the temperature dependence of Xi with 
Ai for v = 3.33 eV A and v — 0, respectively, where 
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FIG. 7: (color online) Vector field plot of the valence-band 
component of d(k) for A2 [(2,2)-component of Eq. (|3T|l ] . 



FIG. 10: (color online) Vector field plot of the conduction- 
band component of d(k) for A4 [(l,l)-component of Eq. I|41[>]. 




FIG. 8: (color online) Vector field plot of the conduction- 
band component of d(k) for A3 [(l,l)-component of Eq. (|3tjp ] . 




FIG. 9: (color online) Vector field plot of the valence-band 
component of d(k) for A3 [(2,2)-component of Eq. 



FIG. 11: (color online) Vector field plot of the valence-band 
component of d(k) for A4 [(2,2)-component of Eq. (|41l) ]. 



v corresponds to the strength of spin-orbit interaction. 
STI with Ai is a full-gapped superconductor, therefore 
resulting Xx> Xy, an d Xz decrease exponentially for T < 
T c both for v = and v = 3.33 eV A. In the case of 
v = 0, all \iS vanish at T = as shown by the dashed 
lines in Fig. [12J On the other hand, in the presence of 
V, all XiS remain to be a finite value at T = (the solid 
line in Fig. [T2J) due to the Van Vleck susceptibility^ 
which is allowed in a multi-band system with spin-orbit 
interaction. Actually, \z at T = is proportional to 
v 2 (See Appendix |A 1| ). It is noted that the value of \z 
is larger than those of \x and Xy m the normal state, 
because of the anisotropy of the energy band. 



D. Anisotropic full-gap A2: Rotation of d- vector 

The d-vector in STI with A2 is parallel to the z-axis 
in the orbital basis. In the absence of spin orbit inter- 
action, the G?-vcctor for A2 in band basis is also parallel 
to the z-axis as shown in Eq. (|37[) . Consequently, only Xz 
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FIG. 12: (color online) Tesmperature dependences of spin susceptibilities \ x (a), \v (b), and \ z (c) of STI with Ai in the 
presence (v = 3.33 eV A, solid line) and absence (v = 0, dashed line) of spin-orbit interaction. The value of X i ls normalized by 
that in the normal state, which is given by X*{T C ) = Xv(T c ) = 0.309x z (T c ) for v = 3.33 eV A and Xx(T c ) = Xv{Tc) = 0.210x z (T c ) 
for v = 0. 
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FIG. 13: (color online) Temperature dependences of spin susceptibilities of STI with A2. 



decreases and vanishes at T = 0, and Xx, Xy are indepen- 
dent of T, as shown in Fig. [T3I In the low temperature, 
Xz is proportional to T since STI with A 2 has a line node 
on the equator for v = 0, as discussed in Sec. HU In the 
presence of v, Xz decreases exponentially for T < T c , as 
denoted by the solid line in Fig. flUt c). In addition to 
this, Xx and Xy slightly decrease [the solid lines in Fig. 
IT51' a)(b)] since the d- vector is rotated so that d x and d y - 
components are induced in the band basis. Xz at T = 
takes a finite value by the following two reasons. First, 
d(k) is not parallel to the z-axis. Second, Van Vleck 
susceptibility arises, as in the case of Ai. 



E. Point node on poles A3: Induced spin-triplet 
pair potential 

For v~0, all x^'s of STI with A3 are independent of 
T, as denoted by the dashed line in Fig. [TU since the 
energy spectrum is gapless (Sec. IIII C[) . On the other 
hand, for v = 3.33 eVA, Xx and Xy decrease down to 
Xl (T = 0)/ X i(T c ) - 0.4 at T = [the solid lines in Fig. 
IT47 a)(b)], while Xz is independent of T [the solid line in 
Fig. fTSTc)]. This behavior is understood from the in- 
duced spin-triplet component d(k) in Eq. due to 



the spin-orbit interaction. The induced (i-vector d(k) is 
parallel to the xy-planc, as shown in Figs. [5]and[5J con- 
sequently the Xz becomes independent of T. Moreover, 
Xx and Xy take a finite value at T = 0, due to the Van 
Vleck susceptibility The similar result has been obtained 
in the bilayer system^ 



F. Point nodes on equator A4: Rotation of d- vector 
and induced spin-singlet pair potential 

For v = 0, due to d(k) \\ x in the band basis, only 
Xx decreases for T < T c and vanishes at T = [the 
dashed line Fig. |15f a)]. In the low temperature, Xx is 
proportional to T since the energy spectrum has a line 
node on the equator (See Sec. [TTJ) . For v = 3.33 eVA, 
Xx decreases with T and is proportional to T 2 in the 
low temperature, except for the residual value at T = 0. 
This residual spin susceptibility comes from the rotated 
d-vector and Van Vleck susceptibility, due to spin-orbit 
interaction. Xz slightly decreases for T < T c since d z (k) 
is present. On the other hand, d y (k) vanishes up to the 
first order of k [Eq. (gT|) and Figs. fTD|[lT]. and thus Xy 
is almost independent of T [the solid line in Fig. fTST b)]. 
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V. DISCUSSION AND SUMMARY 

In this paper, we have calculated the temperature de- 
pendence of the specific heat and the spin susceptibil- 
ity. The temperature dependences of the specific heat 
are similar among three of four possible pair potentials. 
On the other hand, wide varieties of temperature depen- 
dence appear in spin susceptibility depending on the di- 
rection of the applied magnetic field. These results are 
summarized in Table ITTT1 

Finally, we compare the obtained results and the ex- 
perimental ones. From the temperature dependence of 
the specific heat, Ai, A2, and A4 are almost consis- 
tent with the experimental result. On the other hand, 
recent tunneling spectroscopy for the (111) surface of 
Cua;Bi2Se3 has shown a pronounced ZBCPi^ 3 - From the 
theoretical calculation, the ZBCP due to ABSs is gen- 
erated on the (111) surface only for the case with A2 
and A4i 29 i 33 Based on these background, promising pair 
potentials are A2 and A4. In the light of the obtained 
spin susceptibility in this paper, we conclude that it is 
possible to distinguish between A2 and A4 by measuring 
the temperature dependence of in-plane and out-of-planc 
Knight shifts. 



pairing 


specific heat 


Andreev 


spin susceptibility 


potential 


c s 


bound state 


Xx Xy Xz 


Ai = Act 


yes 


no 


/ S v/ 


A2 = ACTyS z 


yes 


yes 


- - v/ 


A 3 = Aa z 


no 


no 


s >/ - 


A 4 = Aa y s x 


yes 


yes 


s - - 



TABLE III: [First column] Possible pairing symmetry of a 
STI. [Second column] Comparison of line shape between our 
results and the experimental one4£ [Third column] Presence 
or absence of ZBCP due to ABSs3>22 [Fourth column] The 
temperature dependence of Xx, Xy: an d Xz- \/ denotes de- 
creasing of Xi with decreasing the temperature. — denotes 
that Xi is almost independent of the temperature. 
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Appendix A: Effects of spin-orbit interaction 

1. Van Vleck susceptibility for Aj 

Here, we derive Van Vleck susceptibility, which gives 
a finite value of the spin susceptibility at T = 0. Wc 
focus on a STI with Ai, based on the Hamiltonian in the 
continuum limit as 

H(k) = [c(k) + m(k)(T x + v z k z a y + vk\\h s (k)a z ]T z 

+ At x . (AI) 

First, we diagonalize spin-part: h s (k) = (k x s) z /k\\, 
where fc|| = |fc|| | = (k 2 + k 2 ) 1 / 2 . The eigenvalue s of h s is 
given by s = ±1. The corresponding eigenvector is given 
by 



(A2) 



with sintpk = fc^/fcy and cosifk = —k y /k\\. 

Next, we diagonalize the normal part: Hos(k) = 
m(k)a x + v z k z a y + svk\\a z . The eigenvalue of Hos(k) 
is given by ar)(k) + c(fc) with 

7?(fc) = yjm 2 (k) + v 2 z k 2 z + v 2 k 2 , (A3) 

and a = ±1 being band index. The corresponding eigen- 
vector \sa) of H os (k) is given by 



l±,±) = 



l±,T) 



cospfe/2 \ 
±e^sm Pk /2j ' 

sinp fc /2 \ 



with 



cosp fc = vk\\/r](k), 

s'mp k = y / m 2 (k) + v 2 k 2 jr\{k), 



cos qk = m(k) / \/m 2 {k) + v z k 2 , 
smq k = v z k z /y/ m 2 (k) + v 2 k 2 . 



(A4) 
(A5) 

(A6) 
(A7) 
(A8) 
(A9) 



In the band basis, the original Hamiltonian is rewritten 
as 

Hsa(k) = [ar)(k) + c(k)}r z + At,. (A10) 
The energy eigenvalue of Hgg-(k) is given by rE^(k) with 



E a {k) = ^[S#)+c(fc)] 2 + A 2 , (All) 

and r = ±1. The corresponding eigenvector \sar) is 
given by 



\s,a,+) = 



s,a, 



cosP fe(f / 2 | 
K smP k */2j ' 

' sinP fe 5/2 \ 
-cosP fc( j/2/ ' 



(A12) 
(A13) 



with cos P k a — [^(fc) + c(fc)]/ 'E„{k). 

For a full-gapped system, spin susceptibility at T = 
is given by 

_2^| ^ \(~s\ Si \s')(~s,a\~s',a')(s,g,+\~s',a',-}\ 2 



E a {k) + E a ,{k) 



(A14) 



For simplicity, we assume that all the g-factors are equal 
to two, and we concentrate on Xz- The matrix element 
in the above expression is estimated as follows. 



(s\s z \s') = 1 - 5 S s>, 
(s,±\ - 5, =f) = ±cosp fc , 



(s, (7, + — s, g , — ) = sin ■ 



Pfecr' Pfccr 



(A15) 
(A16) 

(A17) 



From the above equations, only the Van Vleck term, 
which comes from the off-diagonal term of s ^ s', a =/= a' , 
and t/t', can be nonzero as 



_ 8/i| x - 

c * - iv A, 



l 



v 2 k 2 . Pfe+ _ Pk _ 



N ^ E+(k) + E-(k) T] 2 (k) 2 

A; 



(A18) 



One can verify Xz — > at v — > from the above ex- 
pression. Spin susceptibility at T = in STI with Ai 
stems from Van Vleck component due to the spin-orbit 
coupling v. 



2. Rotation of d-vector for A2 

Here, we derive d-vector for A2 in the band basis. The 
following relation is useful: 

(s,a\a y \s',a') = [S ss > (a z sinp k sinq k - a y cosq k ) 
- (s z ) ss >a x cos p k sinq k 
+ (s x )ss> (5* sinq k -v v sin p k cos q k ) 
+ (s y ) ss >a cos Pk cos qn] ( A1 9) 



The above expression is derived with the use of Eqs. (|A4j) 
and (|A5|) . The pair potential A2 is represented in the 
band basis: 



Aa y s z = A(sinq k o- z - cos q k sinp k a y )s x 
vk\\ 

+ A cos qk—p-Sy, 



(A20) 



where §i is Pauli matrix in the spin-helicity space. Here, 
the relation between s and s is as follows. 

(s\s z \s') = 1 - 8 SS > = (5») g y (A21) 
(s|s x |a') = (s y )ss' sinifk + {s z ) S s' cos^fe, (A22) 
(s\s y \s') = (s z )ss' sinipk - {s y )ss' cosip k , (A23) 
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or equivalcntly, 



Sz, 



s x sin ip k - s y cos (p k 



s x cos (pk + s y sin ipk 



(fc x s) 2 



(A24) 
(A25) 

(A26) 



Consequently, the d-vector for A2 in the band basis is 
obtained as 



d x {k) = Acosqk 
d y {k) = Acosgfc 



vk x 

Vky 

rj(k) ' 



d z {k) = A(smq k a z - cos q k sin p k a y ). 



(A27) 

(A28) 
(A29) 



It is noted that the spin in the above expression is rep- 
resented in the original spin space (s) not in the spin- 
helicity space (s). For STI with A2, spin-orbit interac- 
tion takes a role to rotate d-vector in the band basis. In 
the case of v — 0, due to d(k) \\ z, only Xz decreases for 
T < T c . In the case of v ^ 0, d x (k) and d y {k) (propor- 
tional to v) are present, and Xx and Xy a ls° decreases 
slightly for T <T C . 



which implies that a spin-triplet pair is induced in the 
band basis. Note that d(k) _L z. This is the reason why 
almost only Xz in STI with A3 decreases with tempera- 
ture. 



4. Induced spin-singlet pair and rotation of 
d-vector for A4 



In this subsection, we derive d- vector for A4 = a v s x , 
and show that a spin-singlet pair is induced and the d- 
vector is rotated in the band basis. From the matrix 
elements of a y [Eq. (|A19j) ] and of s x [Eq. (|A22|) ]. the pair 
potential A4 in the band basis is represented as follows, 



Aa y s x = A—— sin q k cr x - A-— cos q k s x 
k 

+ A-^ (sin q k a z - sinp fc cos q k a y )s y 
k 

- A-p-(sinp fe sin q k a z - cos q k a y )s z . (A34) 
fci 



Using Eqs. (|Al4)) . (TAl5|) . and (fS26|, the d- vector for A 4 
in the band basis is obtained as 



3. Induced spin-triplet pair for A3 

In the following, we show that a spin-triplet pair is 
induced for A3 in the band basis. As in Appendix IA 21 
we derive d(k) for A3 = Aa z . By using Eqs. (|A4[) and 
(|A5j) . the matrix elements of a z are obtained as 



(s,a\a z \s,a) = (sa z cosp k + a x sinp fc )~~, . (A30) 
Therefore, a z is expressed in the band basis as 



~ , v(k x Sy kyS X *) _ 

a z = sin p k a x H — a z 

V(k) 



(A31) 



This is derived with the help of Eq. (|A26j) . As a result, 
d (k) and d(k) are given by 

d (k) = Asinp k a x , (A32) 

d(k) = ^a z (-k y ,k x ,0). (A33) 



d (k) = A-^- sin q k a Xl 
r}(k) 



d x (k) = A 



d z (k) 



ky v z k z 



'klmVk) kl \ * 



dy{k)=A kxky 



k 2 
cos q k 



Sin q k 



Vzkz * 

V(k), 



m(k) 



vk x 

-A— — cosq k . 
V{k) 



(A35) 



(A36) 



(A37) 
(A38) 



Therefore, because of spin-orbit interaction, a spin singlet 
pair do(k) is induced in the band basis, and the d- vector 
is rotated so that d y (k) and d z (k) become nonzero. 
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